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An alternative, to the commonly used blast-wave, model describing the freeze-out hypersurface 
is applied to fit the pT-spectra of identified hadrons measured at relativistic heavy-ion collisions at 
^/sNN ~ 62.4, 130 and 200 GeV. Decays of resonances are taken into account completely. It has 
turned out that the fitted freeze-out temperature and baryon number chemical potential depend 
weakly on the centrality of the collision and their values are close to the chemical freeze-out values 
determined from fits to particle yield ratios. 
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I. INTRODUCTION 

During a heavy-ion collision a hot and dense medium is created which eventually evolves into a state of freely 
. streaming particles. The process of hadron decoupling is called freeze-out and two kinds of freeze-out are distinguished 
^ ' Pi 13' (i) chemical freeze- out at Tchem when the hadron abundances become fixed and (ii) thermal (kinetic) freeze- 
, out at Tfci„ when elastic rescattering processes cease and hadrons start to escape freely. And Tchem > Tkin is 
' expected. Values of the statistical parameters at the chemical freeze-out are determined from fits to particle yield 
" ratios, whereas corresponding values at the kinetic freeze-out are fitted to the spectra of hadrons. From here forward, 
statistical parameters at the kinetic freeze-out mean their values fitted to the spectra. Tchem ^ 150 — 170 MeV 
J — \ is estimated at highest heavy-ion reaction energy |3l-[l3|. Additionally, fits done for various centrality classes have 
. revealed that Tchem is almost independent of centralityjl, 0, [sl-fToj . On the contrary, the temperature at the kinetic 
■ freeze-out depends on the centrality and is substantially lower. From the most central to the peripheral bin it changes 
"5 ; as follows: Tkin = 121 - 161 MeV for PHENIX at ,JIWn = 130 GeV Q, Tkin = HI - 147 MeV for PHENIX at 
^ ' = 200 GeV fl5|, Tkin = 89 - 129 MeV for STAR at vs]vjv = 200 GeV @ and Tk^n = HO - 115 MeV for 

O , BRAHMS at ,JIWn = 200 GeV 0. For the PHOBOS data at ./sWn = 62.4 GeV Tk^n = 103, 102, 101 MeV for 
2 ' the central, mid-peripheral and peripheral bin, respectively jl7| . However, the aforementioned estimates oiTkm have 
^ \ been done within the very simplified hydrodynamic model, i.e. the blast-wave model fl8| . 

J> . The main difference between fitting the statistical parameters at the chemical freeze-out and at the thermal freeze- 
out is that the first procedure is independent of a pattern of expansion. Thus, assuming boost invariance and that 
statistical parameters are constant on a freeze-out hypersurface one has 
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where the last equality follows from the factorization of the volume of the system (for more details see [T9l - [2l| ) and 
the density of particle species i reads 



a 

where riP""^°'^'^'^°''-{T, /is) is the thermal density of the ath particle species at the freeze-out, g{i, a) is the final number 
of particles of species i which can be received from all possible simple or sequential decays of particle a and the sum 
is over all species of resonances included in the hadron gas. This means that values of the statistical parameters at 
the chemical freeze-out do not depend on the assumed model of expansion. Therefore fitting measured ratios can not 
distinguish between the blast wave model and the present model. The only sources of valuable information about the 
pattern of expansion are measured spectra. 

In this paper, it will be shown that the behavior of Tkin is model dependent and within a different hypersurface and 
with complete treatment of resonance decays different conclusions about statistical parameters at the kinetic freeze- 
out can be obtained. Namely, the statistical parameters at the kinetic freeze-out are roughly centrality independent 
and their values are close to the corresponding values at the chemical freeze-out. 
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II. FOUNDATIONS OF THE MODEL 

The model applied here is inspired by the single-freeze-out model of Refs. [2^, '23*1, but its main assumption about 
the simultaneous occurrence of chemical and thermal freeze-outs is not postulated. The assumptions of the present 
model are as follows. A noninteracting gas of stable hadrons and resonances at chemical and thermal equilibrium 
appears at the latter stages of a heavy-ion collision. The gas cools and expands, and after reaching the freeze-out 
point it ceases. The conditions for the freeze-out are expressed by values of two independent thermal parameters: 
T and hb- The strangeness chemical potential /zg is determined from the requirement that the overall strangeness 
equals zero. All confirmed resonances up to a mass of 2 GeV from the Particle Data Tables [l^l, together with stable 
hadrons, are constituents of the gas. The freeze-out hypersurface is defined by the equation 



with the additional condition 



const . (3) 



^ Pmax- (4) 



which makes the transverse size of the system finite, t and Pmax constitute two geometric parameters of the model. 
The four- velocity of an element of the freeze-out hypersurface is proportional to its coordinate 



^ ^ 1 (i/j^/jl/^) . (5) 

T T \ t ' t ' t J ' 



The shape of the hypersurface defined by eqs. ©-(IH) and the above form of the flow imply that the model is 
longitudinally boost-invariant and cylindrically symmetric. The last feature is fulfilled strictly only in the case of 
a central collision, but experimental spectra are averaged over the azimuthal angle, which means that they look as 
if they were uniform in this angle. In the view of eq. ([S]), r given by eq. ([3]) is the usual proper time of the fluid 
element. Thus in this model the freeze-out takes place at the same moment of time in the local rest frame of a fluid 
element. In the c.m.s. of the colliding nuclei this means that the farther particles are from the center of the collision 
the later they freeze out. For the 1-1-1 dimensional Bjorken model [25j this is the direct conclusion from the form 
of its solution: all thermodynamical quantities depend only on a local proper time r — \/t'^ — z'^, so if one assumes 
that the freeze-out takes place at a given temperature, T/^.o., it means that it happens at one moment of the proper 
time. The same is true for a spherically symmetric expansion in l-f3 dimensions at sufficiently long times 26]. Note 
that in that case the freeze-out hypersurface has also the shape given by eq. ([3]). The form of the four- velocity, as 
implied by eq. ([5]), determines the three- velocity, v = r/t. This is the so-called scaling solution, which is the exact 
solution for the Bjorken model [25j . In l-|-3 dimensions it was shown that for a^spherically symmetric expansion the 
scaling solution might develop if the sound velocity squared satisfies < 0.2 26]. For boost-invariant, cylindrically 
symmetric systems recent results of ref. ^Jl indicate that the scaling solution can develop in 130 and 200 GeV Au-Au 
collisions within times 7-15 fm. This is roughly the scale of the freeze-out initializing time predicted in this analysis 
for central and mid-central bins (see Table |T|. 

The following parameterization of the hypersurface is chosen: 

i = Tcoshaii coshaj^, r^; = Tsinhaj^ cos(/), = r sinha^ sin(/), = r sinha|| coshaj_, (6) 

where q;|| is the rapidity of the element, a\\ — i&nhr^{rz/t), and a±_ controls the transverse radius r — rsinhaj^. 
The maximum transverse-fiow parameter (or the surface velocity) at the central slice is given by 

pmax _ Pmax _ Pmax / T 

^ + Pmax VI + iPmax/r)^ 



The invariant distribution of the measured particles of species i has the form [22, |23j 

= j p^da,f,{p-u), (8) 
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where da^ is the normal vector on a freeze-out hypersurface, p-u = p^u^ , is the four- velocity of a fluid element and 
fi is the final momentum distribution of the particle in question. The final distribution means here that fi is the sum 
of primordial and simple and sequential decay contributions to the particle distribution. The primordial part of fi is 
given by a Bose- Einstein or a Fermi-Dirac distribution at the freeze-out. A decay contribution is a one-dimensional 
or multidimensional integral of the momentum distribution of a decaying resonance (the exact formulae are obtained 
from the elementary kinematics of a many-body decay or the superposition of such decays, for details see (l9j and 
the Appendix in [28|). The resonance is a constituent of the hadron gas and its distribution is also given by the 
Bose-Einstein (Fermi-Dirac) distribution function. Therefore the final distribution /; depends explicitly on T and fiB- 
With the use of eqs. (O and (O, the invariant distribution ([S]) takes the following form: 

+ 00 Pmax/r 27r 

J da\\ J sinha^c!(sinhQ!x) J d^ p • u fi{p ■ u) , (9) 
-00 



dN, 



d^PT dy 
where 



p ■ u = rriT cosh {a\\ — y) cosh a±_ — pr cos ^ sinh a± . (10) 

The model has four parameters, the two thermal parameters, the temperature T and the baryon number chemical 
potential fiB, and the two geometric parameters, r and Pmax^'t should be stressed that now all parameters of the 
model are fitted simultaneously, opposite to the case of Refs. |22l.[23j where the determination proceeded in two steps. 
First, statistical parameters T and fis were fitted with the use of the experimental ratios of hadron multiplicities at 
midrapidity. Then geometric parameters were determined from fits to the transverse-momentum spectra. Therefore 
the assumption that the chemical freeze-out happens simultaneously with the kinetic freeze-out (the single freeze-out) 
was crucial in that approach. Now all parameters are fitted to the spectra, so the aforementioned assumption is not 
necessary and values of statistical parameters have the meaning of the values at the kinetic freeze-out. 

With the use of eq. ([9]) the measured transverse- momentum spectra of tt^ , , p and p [13, [3 113, HI] can be 
fitted to determine values of the parameters of the model (data points with px > 3 GeV have been excluded) . Fits 
are performed with the help of the method. For the kih measured quantity R^^^ and its theoretical equivalent 
R'j!^{ai, ...,ai), which depends on I parameters ai, ...,a/, the function is defined as 



l,...,a;) = 2^ 2 ' ill-' 



X [a 

fc=i 



where ak is the error of the fcth measurement and n is the total number of data points. The fitted (optimum) values 
of parameters mean the values at which has a minimum. 



III. RESULTS 



The fitted results for Tkin, fJ^B, Pmax and r are collected in Table |I] together with values of the surface velocity 
l^max values of x^/NDF for each centrality class additionally characterized by the number of participants Npart- 
Errors are expressed by one standard deviation (la). Note that except the most peripheral bins of the PHENIX 
measurements all fits are statistically significant. That these new higher values of temperature (higher in comparison 
with the blast- wave fit values) are at minima is shown in Fig. [T] in an example of the temperature dependence of 
for the PHENIX most central bin at ^/snn = 200 GeV. The simulation of the dependence was done within the same 
procedure as in the original single-freeze-out model, with the overall fixed value of hb = 27.5 MeV (the average of 
values from [1, [l^) and the temperature fixed at each point and taken from the range 120-180 MeV with the 1 or 2 
MeV step. Then two geometric parameters were fitted and x^ at their optimum values was put into the figure. 

Coming back to the results of simultaneous fits of all four parameters, some la and 2cr contours are presented 
in Figs. [5] and [3] These contours are done in Tkin and Z?™"^ planes to make them comparable with the blast-wave 
results. But since the parameter space is four dimensional, the n-sigma contour is in fact the 3-dimensional ellipsoid 
embedded in this space, so it can be presented only by means of its projections on some planes. Contours presented in 
Figs. [2] and [3] are such projections onto Tkin and p^"-^ planes chosen at some fixed hb and t, namely at their optimum 
values pB,opt and Topt, Fig. [21 and at pB = PB,opt ± 0.5cr^ and r = Topt ± O.Scr,-, Fig. [31 Note that in Fig. [31black dots 
denote the projections of the points of optimum values of all parameters onto the appropriate plane Tkin - this 
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TABLE I: Values of the statistical and geometric parameters of the model for various centrality bins fitted with the use of 
the RHIC final data for the pr spectra of identified charged hadrons [gI. [T3. [l^. [29l| . All data are at midrapidity, except the 
PHOBOS case (first three rows) where data are at j/ = 0.8 17|- 



Au-Au collision 


Centrality 


N part 




Ms 


pmax 


r 


omax 

P± 


xVndf 


NDF 


case 


[%] 




[MeV] 


[MeV] 


[fm] 


[fm] 








PHOBOS at 


0-15 


294.0 


148.52±1.15 


72.60±3.81 


7.84±0.19 


8.35±0.14 


0.684±0.006 


0.995 


72 


Vsivjv = 62.4 GeV 


15-30 


175.0 


149.64±1.42 


69.47±4.13 


6.23±0.19 


7.20±0.17 


0.654±0.007 


0.43 


72 




30-50 


88.0 


151.29±1.70 


67.15±4.47 


4.60±0.18 


6.05±0.18 


0.605±0.009 


0.20 


71 




0-5 


347.7 


166.74±3.96 


35.06±8.97 


6.31±0.41 


8.08±0.44 


0.616±0.014 


0.53 


78 


PHENIX at 


5-15 


271.3 


161.70±3.21 


43.14±7.77 


6.34±0.35 


7.57±0.34 


0.642±0.012 


0.46 


78 


Vsivjv = 130 GeV 


15-30 


180.2 


162.33±3.29 


38.52±7.75 


5.32±0.29 


6.54±0.29 


0.631±0.012 


0.50 


78 




30-60 


78.5 


162.25±3.46 


31.80±7.87 


3.77±0.23 


4.95±0.23 


0.606±0.015 


0.75 


78 




60-92 


14.3 


159.46±6.85 


37.05±16.09 


1.87±0.27 


3.26±0.27 


0.496±0.043 


1.37 


42 




0-5 


351.4 


150.07±1.34 


24.10±3.66 


9.28±0.21 


9.48±0.19 


0.699±0.004 


0.69 


122 


PHENIX at 


5-10 


299.0 


150.18±1.35 


23.48±3.65 


8.75±0.20 


8.80±0.18 


0.705±0.004 


0.50 


122 


Vsivjv = 200 GeV 


10-15 


253.9 


150.16±1.35 


22.75±3.65 


8.25±0.19 


8.20±0.17 


0.709±0.004 


0.37 


122 




15-20 


215.3 


150.00±1.36 


22.38±3.65 


7.80±0.18 


7.69±0.16 


0.712±0.004 


0.37 


122 




20-30 


166.6 


149.59±1.31 


24.03±3.47 


7.13±0.16 


6.96±0.14 


0.716±0.004 


0.45 


122 




30-40 


114.2 


149.79±1.36 


23.78±3.56 


6.14±0.14 


6.03±0.12 


0.713±0.004 


0.66 


122 




40-50 


74.4 


148.53±1.40 


22.52±3.71 


5.28±0.13 


5.27±0.11 


0.708±0.005 


0.89 


122 




50-60 


45.5 


147.75±1.51 


22.02±4.03 


4.38±0.12 


4.55±0.10 


0.694±0.005 


0.96 


122 




60-70 


25.7 


144.57±1.65 


21.63±4.56 


3.63±0.11 


3.91±0.09 


0.680±0.006 


1.12 


122 




70-80 


13.4 


141.77±1.98 


24.13±5.68 


2.84±0.10 


3.22±0.09 


0.662±0.008 


1.23 


122 




80-92 


6.3 


140.62±2.46 


14.29±7.12 


2.24±0.10 


2.77±0.09 


0.630±0.011 


1.13 


122 




0-5 


352.0 


159.99±1.19 


24.00±2.17 


9.22±0.31 


7.13±0.19 


0.791±0.006 


0.30 


71 


STAR at 


5-10 


299.0 


160.58±1.16 


24.97±2.17 


8.34±0.28 


6.75±0.18 


0.777±0.006 


0.27 


71 


^SNN = 200 GeV 


10-20 


234.0 


161.20±1.14 


22.91±2.15 


7.45±0.24 


6.17±0.16 


0.770±0.006 


0.22 


73 




20-30 


166.0 


162.27±1.12 


23.05±2.17 


6.31±0.20 


5.60±0.14 


0.748±0.007 


0.25 


75 




30-40 


115.0 


161.97±1.08 


20.43±2.17 


5.38±0.17 


5.15±0.12 


0.722±0.007 


0.19 


75 




40-50 


76.0 


162.97±1.08 


21.01±2.21 


4.46±0.14 


4.64±0.11 


0.693±0.008 


0.13 


75 




50-60 


47.0 


163.41±1.07 


18.75±2.25 


3.67±0.12 


4.13±0.10 


0.664±0.008 


0.13 


75 




60-70 


27.0 


162.39±1.06 


16.47±2.31 


2.95±0.10 


3.79±0.09 


0.614±0.010 


0.26 


75 




70-80 


14.0 


163.70±1.15 


15.84±2.50 


2.22±0.09 


3.16±0.08 


0.574±0.012 


0.61 


75 




0-10 


328.0 


150.60±1.39 


23.07±3.51 


9.26±0.25 


8.65±0.21 


0.731±0.004 


0.43 


114 


BRAHMS at 


10-20 


239.0 


151.38±1.48 


26.53±3.72 


8.07±0.23 


7.68±0.19 


0.724±0.005 


0.42 


114 


Vsivjv = 200 GeV 


20-40 


140.0 


149.43±1.54 


25.92±3.98 


7.00±0.21 


6.73±0.17 


0.721±0.005 


0.26 


112 




40-60 


62.0 


148.36±2.02 


26.69±5.21 


5.02±0.20 


5.38±0.17 


0.683±0.008 


0.52 


112 



is the reason why they are outside the contours. In the right panel of Fig. [21 all examples are for most peripheral bins, 
except the PHENIX case at y^ijvw = 130 GeV which is for the second from the most peripheral class. This is because 
errors in the case of the most peripheral bin of the PHENIX measurements at i/F/vjv = 130 GeV are substantially 
greater (see Table |l|, which results in too big extension of la and 2a contours 3 times bigger than in the second 
from the most peripheral bin). This fact together with x^/NDF much higher than 1 does not enable to determine 
firm optimum values of Tkin and in this case. Generally, is flatter in the vicinity of the optimum points for 

PHENIX at ^/sNN = 130 GeV, as one can notice from Figs. [5] and |31 This is also expressed by ^ 2 times greater 
errors of the fitted values in this case. For other experiments la and 2a contours are relatively small. Their sizes do 
not change visibly from most central to mid-central bins and then increase gradually up to the sizes of the contours 
for most peripheral bins (in fact this behavior reflects the behavior of the corresponding errors, cf. Table U]). 

The feeding from weak decays is treated in the same way as experimental groups do, except the STAR case. So, 
for the PHENIX measurements at ^/snn = 200 GeV protons (antiprotons) from A (A) decays are excluded. In 
the PHOBOS case protons (antiprotons) from A and (A and E^) decays are not counted. As far as STAR is 
considered, this Collaboration has corrected its pion spectra for weak decays @. However, STAR also claims there 
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FIG. 1: Simulation of the temperature dependence of x^/NDF for PHENIX at ^sJTn" = 200 GeV and the 0-5% centrality 
class. Solid line is the polynomial 3 best fit. 
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FIG. 2: The contours {la and 2a) in the parameter plane Tkin and /J™"^^ fixed by taking /ib and r at their optimum values. 
In the right panel the PHENIX case at ^ysNN ~ 130 GeV is represented by the second from the most peripheral bin, see text 
for explanations. Dots denotes the optimum values of Tki„ and /J™"^ . 



that when the full feeding from weak decays is taken into account the single-freeze-out model with all the resonances 
can fit the spectra but with higher x^/NDF. Therefore, to check this statement all weak decays are included in the 
STAR data analysis. 

Results for Tkin and /is are also depicted as functions of Npart in Figs. S] and [51 respectively. It is clearly seen that 
both Tkin and b are ahiiost independent of the collision centrality, only for peripheral bins some dependence can be 
observed. Additionally, their values are very close to the values at the chemical freeze-out. Namely, Tchem — 165 — 169 
MeV and = 33 — 38 MeV from the peripheral to most central bin at ^/snn = 130 GeV was found in Ref. [6|, 
Tchem ~ 155 MeV and « 26 MeV independent of the centrality for PHENIX at ^/sWn = 200 GeV in Ref. [s*] and 
Tchem ~ 160 MeV independent of the centrality and /i^ — 15 — 24 MeV from the peripheral to most central bin for 
STAR at ^^sWn = 200 GeV in Refs. d^lll. 

In Figs. [MTOl measured and predicted spectra of identified hadrons are presented for some collisions at RHIC. In 
Fig.[5]the case of PHOBOS at y^sjvjv = 62.4 GeV is dealt with. The way the data and results are depicted is suggested 
by the PHOBOS analysis ^3 (see Fig. 7 there). The lines are predictions of the present model done with the use of 
fitted parameters tabulated in Table IH But fits were done within intermediate ranges ofpT, that is in the ranges where 
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FIG. 3: The contours (la and 2a) in the parameter planes Tkin and /J™"^ fixed by taking fiB and r at ±0.5ct from their 
optimum values, At = t — Topt, A/ns = Hb — HB,opt- AH cases represent the most central classes. Dots denote the projections 
of the points of optimum values of all parameters onto the appropriate plane Tkin - p±"'^. 
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FIG. 4: Centrality dependence of the kinetic freeze-out temperature for the RHIC measurements at ^/snn = 62.4, 130 and 
200 GeV. The lines connect the results and are a guide. 
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FIG. 5: Centrality dependence of the baryon number chemical potential for the RHIC measurements at ^Jsnn = 62.4, 130 and 
200 GeV. The lines connect the results and are a guide. 
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FIG. 6: Invariant yields as a function of pr in Au+Au collisions at y^sjvjv ~ 62.4 GeV. Lines are predictions of the present 
model and symbols are PHOBOS data [17j . 



single charged hadron data exist. Then predictions were made for sums of negatively and positively charged hadrons 
of the same kind but in the whole accessible ranges, namely ranges extended to the very low px where separate data 
on the summed yields exist. As one can see from Fig. [S] predictions of the model agree very well with the data, only 
slight overestimation of protons and antiprotons can be observed but for the very low pt region results agree within 
errors. Predictions of the blast- wave model also agree with the data (cf. Fig. 7 in fTT*!), but opposite to the present 
analysis, they underestimate pion yields and correctly estimate proton and antiproton yields at very low pT- For the 
collisions at y^SNN = 130 GeV, Fig. [71 predictions for pion and kaon spectra agree very well with the PHENIX data, 
whereas for protons and antiprotons the agreement holds up to px ~ 2.5 GeV. Similar results were obtained within 
the blast-wave model (cf. Fig. 19 in JJ.]). In Figs. ISlfTUlresults and data are presented for the PHENIX measurements 
at y'sNN = 200 GeV. Predictions for pions. Fig. [51 agree very well with the data practically in the whole pT range 
for the first five bins from the top, for the next bins the predictions start to miss the data at pt ~ 2.5 GeV and this 
value is decreasing with the centrality to pt ~ 2 GeV for the most peripheral bin. Kaons are reproduced very well in 
the whole range of px and for all centralities, as one can see in Fig. [HI Protons and antiprotons agree very well up to 
Pt ~ 3.3 GeV for first nine bins from the top, for the two last peripheral bins the agreement is lost at pj- sa 2.5 GeV, 
as it can be seen in Fig. [TUl 

A simple test of the model can be performed with the use of the spectra of fl hyperon. This is because D, has 
only the thermal contribution to the invariant distribution, eq. (|9]). Such a test has been done for the blast-wave 
model [lo| and for the single- freeze-out model [l2l |. but in both cases the comparison was done with preliminary 
data. Present results together with the STAR data for il^ + production at ^sj^n ~ 200 GeV are shown in 
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FIG. 7: Invariant yields as a function of pr in Au+Au collisions at ^/snn = 130 GeV. Lines are predictions of the present 
model and symbols are PHENIX data [l^l. For clarity, the data points are scaled vertically for two bins as quoted in the figure. 
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FIG. 8: Invariant yields of tt (left) and ■k~ (right) as a function of pT in Au+Au collisions at -^sjviv = 200 GeV. Lines are 
predictions of the present model and symbols are PHENIX data 29]. For clarity, the data points are scaled vertically for five 
bins as quoted in the figure. 
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Pt [GeV] Pt [GeV] 

FIG. 9: Invariant yields of A"*" (left) and (right) as a function of pr in Au+Au collisions at ^/snn = 200 GeV. Lines are 
predictions of the present model and symbols are PHENIX data [2^. For clarity, the data points are scaled vertically for five 
bins as quoted in the figure. 
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FIG. 10: Invariant yields of protons (left) and anti-protons (right) as a function of pr in Au+Au collisions at y^ijvjv = 200 
GeV. Lines are predictions of the present model and symbols are PHENIX data [2y]. For clarity, the data points are scaled 
vertically for five bins as quoted in the figure. 
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FIG. 11: Transverse momentum distributions of Sl^ + f2 for | y |< 0.75 in Au-Au collisions at y^ijvjv = 200 GeV. Data are 
from 31] (STAR) scaled for clarity, (statistical) errors are of the size of symbols. Lines denote model predictions: solid based 
on fits to the STAR spectra, dashed based on fits to the PHENIX spectra. 



Fig. [TT] Values of parameters for 20 — 40% and 40 — 60% centrality bins explored by STAR in fl measurements 
are the averages of the values from Table |T] for bins which added percent coverage equals 20 — 40% and 40 — 60%, 
respectively. One can see that predictions based on fits to PHENIX spectra agree well with the data. Predictions 
based on fits to STAR spectra agree only qualitatively, they have higher normalization. Also for the — 5% bin the 
slope differs. Blast-wave model predictions for ft^ + 0.^ spectrum for the — 5% bin of the preliminary STAR data 
were done in Ref. [loj . but they do not agree with the data neither in normalization nor in a slope. The probable 
reason for the worse agreement of the STAR data based predictions in the present model, Fig. [11] is that STAR 
spectra of identified stable hadrons are measured in narrower ranges of pr than PHENIX ones, i.e. pr G [0.2,1.2] 
GeV/c approximately for STAR ^ whereas for PHENIX pr € [0.25,2.95] GeV/c (pions), pr £ |0. 55, 1.95] GcV/c 
(kaons) and pt G [0.65,4.25] GeV/c [(anti)protons] [2^. And the STAR measurement of 17^ + 17+ is done within 
the range pr € [1.25,4] GeV/c, practically outside the STAR range of pr of identified stable hadrons but covering 
the great part of PHENIX pT ranges. Also they differ in common ranges of pr, namely STAR spectra are placed 
slightly above the corresponding PHENIX spectra and in the case of pions have different slopes (it has been checked 
carefully for the common — 5% bin after conversion of STAR spectra from tjit — rrii to pt, see Fig. [T2|) . Actually, 
the difference in pion spectra is even greater than what can be seen in Fig. [T^ because the STAR pion spectra are 
corrected for weak decays @. Proton and antiproton spectra are not compared, since the PHENIX subtracted protons 
(antiprotons) from A (A) decays, so by definition the STAR proton (antiproton) spectra are above the PHENIX ones. 
Also that the STAR minimizations have lower x^/NDF than PHENIX ones is probably because it is much easier to 
fit over a narrower range of pr ■ Discrepancies between STAR and PHENIX spectra explain why the optimum values 
of parameters are different for these collaborations. These discrepancies could be caused by the different detectors 
used by the aforementioned Collaborations. The PHENIX uses a spectrometer, which consists of drift chamber (DC), 
pad chamber (PC) and time-of- flight (TOF) 29], whereas the STAR uses time projection chamber (TPC) Q. The 
TPC is better in measurements of resonances but the PHENIX spectrometer measures spectra of identified particles 
more precisely [1, HI]- This might also help to understand why f7~ +0"*" measured spectra (an example of the STAR 
data on resonance production) are well reproduced in the model with the use of the parameters fitted to the PHENIX 
data on identified stable hadrons (in addition to the earlier argument based on the compatibility oi px ranges). 



IV. DISCUSSION AND CONCLUSIONS 



Before the final conclusion one result of ref. [3] should be commented. In Q a hydrodynamical model supplemented 
with the dynamical freeze-out criterium (33ri35i] 



'^scatt '^exp 



(12) 



11 




FIG. 12: Comparison of tt^, tt", A"*" and K~ spectra measured by STAR and PHENIX for the — 5% centrality bin at 
y^SNN ~ 200 GeV. All STAR data points are depicted, whereas PHENIX data ranges are cut from the right side (in fact they 
extend to pr = 2.95 GeV for pious and pr = 1.95 GeV for kaons). 



where 1/Tscatt is the local scattering rate and I/t^xp is the local expansion rate and the parameter ^ ~ 1, was 
successfully used to explain the centrality dependence and the temperature range of T^in determined from the blast- 
wave fits to the STAR data at ^Jsnn = 200 GeV To obtain the T^in range comparable with the STAR results, 
Tkin — 90 — 130 MeV, the appropriate adjustment of the parameter ^ was done, ^ — 0.295. Then, the range 
Tkin — 105 — 135 MeV was obtained in [2!]. According to the author's knowledge the parameter ^ appeared the first 
time in Q just to fix the range of temperature as was explained above. 

Generally, the dynamical criterium for freeze-out was formulated as the moment when the local velocity of the 
rarefaction due to the expansion overcomes the local thermal velocity [1^. Since the freeze-out means the end of 
collisions the criterium can be stated more precisely as the moment at a space point when the average local proper 
time between subsequent collisions Tscatt overcomes the characteristic time for expansion of the system given by 
Texp = V/V, where V denotes the derivative of the volume of the system with respect to the local proper time. The 
expansion characteristic time is the time which the evolution of the system from very small volume (almost zero) to 
the given volume V would take with the constant speed equal to the local velocity V . For a given hadron species 
the scattering time can be estimated as Tgcatt = ^/vth = ^/ PtotO'totVth, where A is the mean free path, Ptnt is the 
total density of all particles, atot is the average cross section and Vth is the average thermal velocity [33, [S^- The 
characteristic expansion time can be expressed as Texp = = V/V [36j . 

Hydrodynamics can be applied for Tgcatt ^ Texp but when Tscatt ^ Texp certainly a system is not in statistical and 
thermal equilibrium 37[ . However the determination of the strict border between a flow and a free stream of particles 
is somehow arbitrary, it is commonly accepted that on this border Tscatt ~ Texp- Both times increase when the system 
cools down, but Tscatt is growing steeper, so it cuts Texp at some temperature (it represents the case ^ = 1; for some 
examples of simulations of Tscatt and Texp see ref. [s^ ) . In terms of rates t^^;^^ j is falling steeper and cuts r^^^ at some 
temperature. For the expansion rate multiplied by ^, if ^ < 1 then fr^^ cuts rZL^ at some lower temperature, if 
^ > 1 it happens at some higher temperature. Therefore, the choice ^ = 0.295 as in [3] means decreasing the freeze-out 
temperature so as it falls into the blast-wave range. Thus it is very likely that increasing ^ but still keeping ^ ^ 1 
one could obtain the freeze-out temperature in the range Tkin — 140 — 165 MeV as in the present model. Another 
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point is the centrality dependence of Tf^m- The change of T^m with the centrahty, of the order of 25% between the 
most central and the most peripheral bins, was obtained from eq. with the same expression for r"^^^ taken for 
all centrality classes 0. Such assumption might not be reasonable, since e.g. in Ref. [s^ variation of the temperature 
dependence of Tscatt with specific entropy S/A = s/ub (entropy density per baryon number density) was shown, see 
Fig. 7 there. The specific entropy, as an initial condition, could be different for each centrality class, so the different 
expression for t"^^^ should be put into eq. in each case. Moreover, the scattering rate was calculated for pions 
only in a pion-kaon mixture [37|. The full scattering rate should include interactions between all constituents of the 
gas. Of course, this is an open question how such realistic scattering rate could influence the centrality dependence 
of Tkin determined from eq. (HU, but that this might flatten the dependence can not be excluded. And the last 
remark, Tkin points described as the STAR data in Fig. 3 of Ref. [1] are not the measured quantities. They are 
blast- wave- fit results. Thus Fig. 3 in Q demonstrates that T^m obtained within the model where the freeze-out takes 
place on a hypersurface of constant temperature is consistent with the average T^in for the hypersurfacc on where the 
temperature is not constant. 

In summary, the alternative hydrodynamical model has been proposed to describe hadronic px spectra measured at 
relativistic heavy-ion collisions. So far, conclusions about chemical and thermal (kinetic) freeze-outs have been drawn 
from the blast-wave parametrization |18l] of the final stage of the collision . It has turned out that those conclusions are 
not definite and depend strongly on the applied model. In the present model the temperature and the baryon number 
chemical potential fitted to the spectra are almost independent of the centrality of the collision and their values are 
very close to the values at the chemical freeze-out, what is opposite to the conclusions drawn from the blast-wave 
model analysis |2i]. Such behavior justifies the ad hoc assumption about one freeze-out postulated in Refs. [2^ [23|. 
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